Suppression of Charge Equilibration leading to the Synthesis of Exotic Nuclei 
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Charge equilibration between two colliding nuclei can take place in the early stage of heavy-ion 
collisions. A basic mechanism of charge equilibration is presented in terms of the extension of single- 
particle motion from one nucleus to the other, from which the upper energy-limit of the bombarding 
energy is introduced for significant charge equilibration at the early stage of the collision. The for- 
mula for this limit is presented, and is compared to various experimental data. It is examined also by 
comparison to three-dimensional time-dependent density functional calculations. The suppression of 
charge equilibration, which appears in collisions at the energies beyond the upper energy-limit, gives 
rise to remarkable effects on the synthesis of exotic nuclei with extreme proton-neutron asymmetry. 

PACS numbers: 21.10.Ft, 25.70.-z, 25.70.Hi 



Charge equilibration is a rapid process during the early 
stage of heavy-ion collisions with a time scale of 10 -22 sec 
(for a review, see [l|). Despite many theoretical at- 
tempts, its mechanism has remained an open problem, 
where only the relation with the isovector giant dipole 
resonance (iv-GDR) was discussed in certain cases (for 
example, see The charge equilibration is quite 

important, because it naturally prevents the synthesis of 
exotic nuclei with extreme proton-neutron asymmetry. 
As there are and will be the 3rd generation Rl-beam fa- 
cilities, it is an urgent question whether such synthesis 
can be enhanced with higher beam energies or not. 

In this Letter, we first point out that there is an upper 
limit of the bombarding energy for the fast and significant 
charge equilibration in the initial stage of the collision. 
Microscopic three-dimensional (time-dependent) density 
functional calculations are systematically performed for 
the purpose of examination. We actually employ Skyrme 
time-dependent Hartree-Fock (TDHF) theory, which is 
a rather unique presently feasible method for the treat- 
ment of non-perturbative processes such as multi-nucleon 
transfers in a realistic framework. 

The suppression of charge equilibration brings about a 
favorable situation for the synthesis of exotic nuclei. The 
evaluation of the upper limit thus can have crucial sig- 
nificance for experiments on nuclei with extreme proton- 
neutron asymmetry. We shall consider the collision of a 
target nucleus with mass number A\, neutron (proton) 
number N\ (Z\), with a projectile nucleus with A2, N%, 
and Z 2 . The total mass, neutron, and proton numbers 
are denoted by A, N, Z, respectively. We look at this 
problem from an intuitive and basic viewpoint. We begin 
with a picture that the charge equilibration takes place 
as wave functions of nucleons propagate from their orig- 
inal nucleus to the other nucleus in the initial stage of 
the collision. Namely, the regime of individual single- 
particle motion spreads out following the lowering of po- 
tential barrier between the two nuclei after the touching. 



Because this spreading occurs as a consequence of un- 
blocked single-particle motion, the process can be very 
fast; a particle travels into the other side within ~ 10~ 22 
sec at a quarter of light velocity, which roughly corre- 
sponds to the Fermi energy of normal nuclear matter. 
On the other hand, this propagation can be prohibited 
by the Pauli principle, while the neutron (or proton) ex- 
cess weakens this effect. Another important factor is the 
difference of velocities of the two nuclei. Because this 
spreading needs a certain time, it does not lead to charge 
equilibration in the initial stage if the relative velocity of 
the colliding nuclei is too large. In the charge equilibra- 
tion, protons and neutrons from both nuclei, particularly 
those near the Fermi levels, are mixed within a certain 
time after the touching. 

We can introduce an ansatz that, in order for the fast 
charge equilibration to occur, the relative velocity v r of 
the two nuclei at the collision must be below the velocities 
corresponding to the proton or neutron Fermi momenta 
of both nuclei. By denoting the minimum of these four 
velocities as Up 3111 , the present effect can be summarized 
by the statement that the upper limit of the bombarding 
energy for charge equilibration is determined by v r = 
w^ 11111 . The upper limit of the energy in the laboratory 
frame for charge equilibration is expressed as the sum 
of the kinetic energy for velocity v^ 1111 and the Coulomb 
energy at touching: 
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FIG. 1: Yield distribution of final fragments as a function of their N/Z ratios for the collisions of 20S Pb + 132 Sn (N/Z ratio 
is discretized by 0.05) based on TDHF calculation (SLy4d). Four cases with different E C m/A values are presented. Columns 
corresponding to the equilibrium value of 7V/Z=1.58 are colored in blue, and connected by red lines. For reference, the N/Z 
ratios for 208 Pb and 132 Sn are 1.54 and 1.64, respectively. 



TABLE I: E cm /A values [MeV] of the upper limit of charge 
equilibration obtained by TDHF calculations with SLy4d and 
SkM* parameters, compared to those by the proposed for- 
mula, eq. |T]). For reference, the values obtained by the Fermi 
gas model with standard parameters are shown. 





Collision 




TDHF 


TDHF 


Eq. © 


Fermi 








(SLy4d) 


(SkM*) 




gas 


(i) 


' MH Pb + 


238-q 


6.5±0.5 


6.5±0.5 


6.91 


9.46 


(ii) 


208p b + 


132 Xe 


6.5±0.5 


6.5±0.5 


6.50 


9.03 


(iii) 


208p b + 


132 Sn 


6.5±0.5 


6.5±0.5 


6.36 


9.03 


(iv) 


20Sp b + 


40 Ca 


3.5±0.5 


3.5±0.5 


3.66 


5.14 


(v) 


208p b + 


24 Mg 


2.5±0.5 


2.5±0.5 


2.36 


3.52 


(vi) 


208p b + 


24q 


2.5±0.5 


2.5±0.5 


2.18 


3.52 


(vii) 


208p b + 


16 o 


1.5±0.5 


1.5±0.5 


1.75 


2.50 


(viii) 


208p b + 


4 He 


< 1.0 


< 1.0 


0.48 


0.70 


(ix) 


24 Mg + 


24q 


5.5±1.0 


5.5±1.0 


5.99 


9.5 



unit, the vacuum permittivity, and the usual nuclear ra- 
dius parameter (1.2 fm), respectively. Here we express 
the minimum velocity by the corresponding minimum of 
the proton or neutron density of the two nuclei using the 
formula w F min = h(3ir 2 p min ) 1/3 /m. In Eq. ©, i = 1,2 
distinguishes the two initial nuclei, and e and 5 (func- 
tions of Ai) are introduced based on the droplet model 
0, [l(J to take into account the effect of neutron and pro- 
ton skins. 

In order to confirm the validity of this picture, we 
perform three-dimensional TDHF calculations systemat- 
ically. Many TDHF events were obtained with different 
values of impact parameter being incremented by 0.25 
fm. These events are summed over impact parameters 
with weights of geometric cross section, and elastic cases 
are discarded. By collecting TDHF events, we sort them 
in terms of N/Z ratio with N/Z being discretized into 
bins of width 0.05. Figure Q] shows the yield distribu- 
tion of final fragments for 208 Pb + 132 Sn reaction. Go- 
ing from low to high E cm (total kinetic energy in the 
center-of-mass frame) of the collision, the peak energy 
of the yield distribution as a function of N/Z is almost 



constant at the beginning, and is shifted later with the 
lowering of peak height. A clear decrease of the yield 
of charge-equilibrated fragments for E cm /A > 7.0 MeV 
is noticed. On the other hand, very neutron-rich nuclei 
with N/Z ~ 2.0 simultaneously start to be produced. By 
taking E\ as the energy at which the peak height starts 
to be lowered (6.0 MeV in this case), and E2 as the en- 
ergy at which the yield of the equilibrated N/Z becomes 
about 50 % (7.0 MeV in this case), the upper energy-limit 
is defined in TDHF as E C e = (Ei + E 2 )/2 (6.5 MeV in 
this case). The uncertainty from the energy bin value is 
(E 2 -E 1 )/2 (0.5 MeV in this case). Such TDHF results 
are summarized in the third and forth columns of Table 
HI TDHF calculations with two different parameter-sets 
result in the same upper energy-limit. Note that a large 
value is obtained for 24 Mg + 24 0, implying that a sim- 
ple mass dependence (small values for reactions between 
light nuclei, and vice versa) cannot explain this, while 
the mass asymmetry plays a certain role. 

Let us have a comparison between the present upper- 
limit and the TDHF results. Figure [2] depicts how the 
upper limit given by Eq. (Q]) changes as functions of 
N/Z and Ai/A 2 (or A 2 /Ai). The upper energy-limit 
comes down significantly low for higher mass asymme- 
try, while it depends only weakly on the total mass. Al- 
though charge equilibration can compete with Coulomb 
excitation particularly in collisions involving nuclei with 
large Z, no evidence of a major change due to large Z 
is found in the results of Eq. ([I]) or those of TDHF cal- 
culations. However, non- negligible decrease of the upper 
energy-limit due to the proton-neutron asymmetry is no- 
ticed. The corresponding values obtained by the formula 
in Eq. {T]) is shown in the fifth column of Table HI For 
reference, the values obtained by the simple Fermi gas 
model (kp = 1.36 fm -1 for both protons and neutrons) 
are shown in its sixth column. By comparing the re- 
sults of Eq. ([1]) with the TDHF results shown in the 
third and forth columns, the agreement is remarkable, 
including a high value in the last row. Comparison be- 
tween the results of Eq. JT]) (blue circles) and the TDHF 
results (blue bars) is also made in Fig. [5] Consequently, 
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FIG. 2: N/Z dependence of the upper energy-limit of charge 
equilibration based on Eq. ([TJ) in the center-of-mass frame, 
where A\ > A2 is assumed without loss of generality. Values 
with different total masses A = 100, 200, 300, 400, and 500 
are plotted for each A1/A2, which correspond to black lines 
from bottom to top in each group (values are too crowded to 
distinguish the total mass difference, for the cases of A1/A2 = 
10 and 100), and lines are drawn to guide eyes. Each TDHF 
calculation is shown as a blue bar, where the central points 
correspond to the value obtained by Eq. ([T]) for each reaction, 
and Roman numbers distinguish reactions shown in Table [I] 



the upper-limit of charge equilibration depends largely 
on the Fermi energy, and the proton-neutron asymmetry 
can contribute to the shift of the upper-limit, where the 
dependence of the upper-limit on the mass-asymmetry is 
remarkable. 

Let us move on to the comparison to the experiments. 
It can be seen that existing experimental data agrees with 
the present upper-limit formula. For instance, the fol- 
lowing experiments show charge equilibration: 40 Ar + 
58 Ni at E lab /A=7.0 MeVEH, and 56 Fe + 165 Ho ( 209 Bi) 
at Ei ab /A= 8.3 MeV [i~2f~^ The following experiment 
shows the disappearance of char ge e quilibration: 112 Sn 
+ 124 Sn at E lab /A=50 MeV [3, [13 . Recently, experi- 
ments: 124 < 112 Sn + 124 < 112 Sn at E lah /A= 35 and 50 MeV 
resp ectively was performed at Michigan State University 
[161 ]. One remarkable result here is that the final frag- 
ments are not so close to charge equilibrium even in lower 
energy collisions with Ei ab /A=35 MeV. The disappear- 
ance of charge equilibration for this experiment has not 
been explained. Because the upper limit of charge equi- 
libration is calculated to be Ei ab / A=27 .6 MeV from Eq. 
([1]), this question can be understood now. 

The suppression of charge equilibration contribute nat- 
urally to the production of exotic fragment; more ex- 
otic nuclei far from the equilibrated N/Z ratio arc to be 
synthesized above the present upper-limit. Because the 
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FIG. 3: Real-time dynamics of charge distribution for Pb 
+ 40 Ca (SLy4d) is shown for a fixed impact parameter 7.5 
fm. The upper and the lower panels show cases with E cm /A 
= 3.0 MeV and 4.0 MeV, respectively. For both cases, 208 Pb 
is coming from the left, and 40 Ca from the right. The colored 
regions show the distribution of charge (proton-rich part), 
and the contours being incremented by 0.02 fm -3 show the 
density. Domains with the proton-to-neutron density ratio 
greater than 0.78 and 0.72 are colored in red and yellow, re- 
spectively. 



bombarding energies of the 3rd generation Rl-beam facil- 
ities are sufficiently high to exceed the upper-limit, the 
present upper-limit ensures more production of further 
exotic isotopes by the latest and the future Rl-beam fa- 
cilities. In fact, in the experiment of Ref. [l5j . the yield 
of exotic fragments was increased simply by putting the 
beam energy higher than the present upper-limit. It is 
also of much interest to explore the novel possibility of the 
synthesis of exotic nuclei by collisions including only /3- 
stable nuclei. Such possibilities have not attracted much 
attention, but may emerge with various feasibilities of 
the production of exotic spieces, if the experiment is set 
for energies beyond the present upper-limit. As an ex- 
ample, Fig. [3] shows the real-time dynamics of 208 Pb + 
40 Ca around the energy where the production of exotic 
nuclei starts; the upper energy limit is 3.66 MeV (See 
Table HJ), fusion appears at E cm /A = 3.0 MeV (upper 
panels), and break-up is seen at 4.0 MeV (lower panels). 
/3-unstable fragments are emitted only for the bombard- 
ing energy higher than the upper-limit; for instance, a 
small fragment in the lower-right panel of Fig. [3] is 62 Mn 
(numbers of nucleon are rounded to be integer), for which 
the stable isotope is 55 Mn. 

We now point out some interesting details of charge- 
equilibration dynamics for collisions below the present 
upper energy limit. We shall first investigate it with a 
focus on the iv-GDR in a collision between light nuclei. 
Figure [4] (a) shows the time-evolution of charge distri- 
bution for 24 + 24 Mg. We see the appearance of the 
iv-GDR and the oscillation of well-localized charge (4.5, 
6.0 xlO -22 sec). The charge equilibration is synchro- 
nized with the iv-GDR. This is, however, seen only for 
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FIG. 4: Real-time dynamics of charge distribution for 24 Mg 
+ 24 and 24 Mg + 208 Pb (SLy4d) are shown for E cm /A = 4.5 
MeV and 2.0 MeV with the impact parameters 5.0 fm and 7.5 
fm, respectively. For both cases, 24 Mg is incoming from the 
left, and the states are evolving into fusion. The description 
manner is the same as Fig. O where parts with the proton to 
neutron density ratio larger than 1.00 (0.75) and 0.80 (0.70) 
are colored in red and yellow in case of 24 Mg + 24 ( 24 Mg + 
208 Pb), respectively. 



the initial stage. Energy-dependent equilibration mecha- 
nisms of a similar kind have been reported in condensed 
matter physics (see fl9l - l2l"l | ) . where the equilibration is 
fast for lower energies (temperatures) , and quite slow for 
higher energies. Such a similarity between completely 
different physical systems is interesting. 

In this Letter the mechanism of charge equilibration, 
in which nucleons with the Fermi velocity play a primary 
role, has been presented, and its validity is examined by 
comparison to virtually all existing relevant experimental 
data. This concept has been further analyzed in terms of 
a three-dimensional time-dependent density function for- 
malism. The upper energy-limit of charge equilibration 
has a crucial impact on the nuclear synthesis of exotic nu- 
clei. Properties presented in this Letter will give a sound 
motivation for the production of further exotic isotopes 
by the latest and the future Rl-beam facilities. 

This work was supported in part by EMMI project 
of the Hclmholtz Alliance, by JSPS Core-to-Core pro- 
gram EFES and by grant-in-aid for Scientific Research 
(A) 20244022. Y.I. thanks for the JSPS fellowship, and 
express his gratitude to Dr. C. Simcnel. T.O. thanks 
Prof. P. Bonche for long-term valuable collaboration. 



collisions between light nuclei. Similar results are ob- 
tained in [H-Q . Next, we move on to collisions involving 
a heavier nucleus; e.g. 24 Mg + 208 Pb, as shown in Fig. I 
(b). Charge-equilibration dynamics is not similar to the 
motion of the iv-GDR in this case, because no significant 
oscillations of the localized charge appear. Instead a radi- 
ally layered structure of the composite nucleus is formed 
(7.5 xlO -22 sec), in which a relatively neutron-rich core 
appears in the center, a proton-rich layer surrounds it, 
and a neutron-rich skin is in the surface. This seems to 
be due to the Coulomb repulsion, hence the radial distri- 
bution of charge is formed. Here one can find analogous 
situation of isovector monopole excitation. It means that 
the radial flow plays a prominent role in such low energy 
collisions. Similar layered structures are obtained in the 
TDHF calculations (with both Skyrme parameter-sets) 
listed in Table U except for 24 Mg + 24 0. Consequently, 
charge equilibration, which is dependent on the mass, has 
been clarified to be related not only with the isovector gi- 
ant dipole excitations, but also with the isovector giant 
monopole excitations. 

Finally, let us comment on the diffusion-type mech- 
anism towards charge equilibrium. The mechanism of 
charge equilibration discussed in this Letter is valid at 
lower energies (lower than the upper energy-limit). At 
much higher energies, collisions between nucleons and 
the diffusion contribute mostly to charge equilibration 
(for example, see [l7, 18|). As the diffusion is a slow pro- 
cess, it is unlikely to attain the charge equilibrium within 
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